


On the role of the final state interactions in rare 

B-decays 

A.B. Kaidalov, M.I. Vysotsky 
bo"' Institute of Theoretical and Experimental Physics, Moscow 

o ■ 

^ : July 14, 2008 

Abstract 

The effects of final state interactions (FSI) in liadronic B-decays are investigated. Tlie 
r-| ! model for FSI, based on Regge phenomenology of high-energy hadronic interactions is 

Ph| proposed. It is shown that this model explains the pattern of phases in matrix elements of 

C^' B —>■ TTTT and B ^ pp decays. These phases play an important role for CP-violation in 

B-decays. The most precise determination of the unitarity triangle angle a from B^ — > pn 
decays is performed. The relation between CP-asymmetries in i? — > Kit decays is 
discussed. It is emphasized that the large distance FSI can explain the structure of 
^ I polarizations of the vector mesons in B-decays and other puzzles like a very large 

^ I branching ratio of the B-decay to ScAc. 
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1 Introduction 

r> I In this paper we give a review of some unusual properties of the matrix elements in the 

c^ I hadronic B-decays. It is based on papers [1, 2], where B -^ rni, B ^ pp decays were discussed 

and it contains some new material on B ^ pn, B ^ Ktt decays and polarization of vector 
mesons in B-decays. The detailed information on B-decays, obtained in the experiments 
at B-factories [3], provides a testing ground for theoretical models. The investigation of 
rare B-decays and CP violation in these decays provides not only the information on CKM 
matrix, but also on QCD dynamics both at small and large distances. 

One of the most interesting and unsolved problems in B-decays is the role of FSI. In 
this paper we shall demonstrate that FSI play an important role in the hadronic B-decays 
and enable to explain some puzzles observed in rare B-decays. In particular it will be 
demonstrated that the phases due to strong interactions are substantial in some hadronic 
B-decays. These phases are important for understanding the pattern of CP-violation in 
rare B-decays. The model for calculation of FSI will be formulated and compared to the 
data on B ^ TiTT and B -^ pp decays. The model is based on Regge-picture for high- 
energy binary amplitudes and enables to explain a pattern of helicity amplitudes in some 
B-decays to vector mesons. The large distance interactions provide a simple explanation of 
the anomalously large branching ratio of the B-decay to ScAc. The CP-violation asymmetries 



will be discussed and the most accurate determination of the unitarity triangle angle a will 
be presented. 



2 B ^ TiTi / B -^ pp puzzle 



The probabilities of three B ^ im and three B ^ pp decays are measured now with a 
good accuracy and presented in Table I. There is a large difference between the ratios of the 
charged averaged B^ decay probabilities to the charged and neutral mesons: 
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It was demonstrated in refs.[l, 2] that this difference is related to the difference of phases 
due to strong interactions for matrix elements oi B ^ titi and B -^ pp-decays. The matrix 
elements of these decays can be expressed in terms of amplitudes with isospin zero and two. 
To take into account the differences in CKM phases for tree and penguin contributions we 
separate the amplitude with 1=0 into the corresponding parts Aq and P: 
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where Vik are the elements of CKM matrix, 7 and j3 are the unitarity triangle angles and we 
factor out the product m^/^/+(0) which appears when the decay amplitudes are calculated 
in the factorization approximation. 

The charge conjugate amplitudes are obtained by the same formulas with substitution 

/5, 7 ->-/?, -7. 

The CP asymmetries are given by [4]: 
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where tttt is vr+vr" or tt'^tt". 

The analogous formulas take place for pp final states where the longitudinal polarizations 
of p-mesons are dominant. 

The values of P can be determined using d ^^ s interchange symmetry from decays B^ — > 
K^*p^ and B.^ -^ K^n^ [5] and turn out to be rather small compared to tree contributions. 
Note, however, that P determines the magnitudes of the direct CP violation in hadronic 
decays. 



If we neglect the penguin contribution, then the difference of phases is expressed in terms 
of the branching ratios as follows 
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Using the experimental information on the branching ratios of i? ^ Trvr-decays [3] we obtain 
\5l-8l\ = 48^. 

The penguin contributions to Bik do not interfere with the tree ones because CKM angle 
Q; = 7r — /5 — 7is almost equal to 7r/2. Taking into account P-term we get: 

\5l-5l\ = 37° ±10° . (6) 

This agrees with the result of the analysis in ref . [6] : 

5- -5- = 40° ±7° . (7) 

Thus the difference of the phases of the matrix elements with 1=0 and 1=2 is not small in 
sharp contrast with the factorization approximation often used for estimates of heavy meson 
decays. 

For B -^ pp-decays we obtain in the analogous way: 

K-6^,\ = irtr,^ . (8) 

This phase difference is smaller than for pions and is consistent with zero. 

The fact that the phases due to FSI are in general not small for heavy quark decays is 
confirmed by the other D and B-decays. The data on D ^ vr+vr", D -^ 7r°7r° and D"^ -^ tt^tt*^ 
branching ratios lead to [7]: 

(5f -5o^ = 5d = ±(86°±4°) . (9) 

The last example is -B ^ Dti decays. Dti pair produced in i?-decays can have / = 1/2 or 3/2. 
From the measurement of the probabilities of B^ -^ D^7r~, B^ -^ D^n^ and B^ -^ D^n^ 
decays in paper [8] the FSI phase difference of these two amplitudes was determined: 

(5,3. = 29° ±4° . (10) 

Thus the experimental data indicate that the phases due to FSI are not small for heavy 
meson decays. 

3 Calculation of the FSI phases of B ^ tttt and B -^ pp 
decay amplitudes 

Let us remind that for K ^ -k-k decays there are no inelastic channels, Migdal- Watson (MW) 
theorem is applicable and strong interaction phases of S'-matrix elements oi K ^ (Stt)/ 
decays are equal to the phases of the corresponding tttt -^ nn scattering amplitudes at 
E = rriK- 

For B-mesons there are many opened inelastic channels and MW theorem is not directly 
applicable. Serious arguments that strong phases should disappear in the Mq -^ oo limit 



were given by J.D. Bjorken [9]. He emphasized the fact that the characteristic configurations 
of the hght quarks produced in the decay have small size ~ 1/Mq and FSI interaction cross 
sections should decrease as 1/Mq. Similar arguments were applied in the analysis of heavy 
quark decays in the QCD perturbation theory [10]. These arguments can be applied to 
the total hadronic decay rates. For individual decay channels (like B -^ tttt) which are 
suppressed in the limit Mq -^ oo the situation is more delicate. However, even in these 
situations the arguments of Bjorken that due to large formation times the final particles are 
formed and can interact only at large distances from the point of the decay seem relevant. 

On the other hand, the formal analysis of different classes of Feynman diagrams, in- 
cluding soft rescatterings [11, 12], show that the diagrams with pomeron exchange in the 
FSI-amplitudes do not decrease as Mq increases. The same conclusions follow from the 
applications of generalizations of MW-theorem [13, 14]. 

In the process of the analysis of FSI in heavy meson decays it is important to understand 
the structure of the intermediate multiparticle states. It was shown in ref.[2] that the bulk 
of multiparticle states produced in heavy meson decays has a small probability to transform 
into two-meson final state and only quasi two-particle intermediate states XY with the 
masses Mjr^y) < Mb^qcd <^ ^1 can be effectively transformed into the final two-meson 
state. In refs. [1, 2] in calculation of FSI effects for B ^ tttt and B ^ pp decays only two 
particle intermediate states with positive G-parity to which 5-mesons have relatively large 
decay probabilities were considered. Alongside with tttt and pp there is only one such state: 

We shall use Feynman diagrams approach to calculate FSI phases from the diagram with 
the low mass intermediate states X and Y. Integrating over loop momenta d'^k one can 
transform the integral over ko and k^ into the integral over the invariant masses of clusters 
of intermediate particles X and Y: 
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and deform integration contours in such a way that only the low mass intermediate states 
contributions are taken into account while the contribution of heavy states being small is 
neglected. In this way we get: 

Mi^ = MPy\6^xS^y + ^Ti=l^J , (12) 

where Mxy are the decay matrix elements without FSI interactions and TxY^_^nn is the J = 
partial wave amplitude of the process XY -^ tttt (T"' = {S"' — l)/(2i)) which originates 
from the integral over d'^k±. 

For real T Eq.(12) coincides with the application of the unitarity condition for the calcu- 
lation of the imaginary part of M while for the imaginary T the corrections to the real part 
of M are generated. 

This approach is analogous to the FSI calculations performed in paper [15]. In [15] 
2—^2 scattering amplitudes were considered to be due to elementary particle exchanges 
in the t-channel. For vector particle exchanges s-channel partial wave amplitudes behave 
as s'^~^ ~ s° and thus do not decrease with energy (decaying meson mass). However it 
is well known that the correct behavior is given by Regge theory: s"»(o)-i. For p-exchange 
ap{0) ~ 1/2 and the amplitude decrease with energy as 1/^/s. This effect is very spectacular 
for B -^ DD -^ TTTT chain with D*{D2) exchange in t-channel: aD*{0) ~ —1 and reggeized 
D* meson exchange is damped as s~^ ^ 10~^ in comparison with the elementary D* exchange 



(see for example [14]). For 7r-exchange, which gives a dominant contribution to pp -^ tttt 
transition (see below), in the small t region the pion is close to mass shell and its reggeization 
is not important. 

Note that the pomeron contribution does not decrease for Mq -^ oo, however it does 
not contribute to the difference of phases 6q — 62 which we are interested in. So this phase 
difference is determined by the secondary exchanges (p, tt) and it decreases at least as 1/Mq 
for large Mq in accordance with Bjorken arguments. For phases 6q and 62 separately the 
pomeron contribution does not cancel in general. If Bjorken arguments are valid for these 
quantities it can happen only under exact cancellation of different diffractively produced 
intermediate states and it does not happen in the model of refs. [1, 2]. 

Let us calculate the imaginary parts oi B ^ -kit decay amplitudes which originate from 
B ^ pp —^ TCTT chain : 

/// rns 
—M{pp^TnT)M*{B^pp) . (13) 
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In the amplitude pp -^ tttt of pp intermediate state the exchange by pion trajectory in 
the t-channel dominates. It was already stressed that p-mesons produced in 5-decays are 
almost entirely longitudinally polarized. That is why it is necessary to take into account only 
longitudinal polarizations for the intermediate p-mesons. The amplitude of p^p^ -^ tt^tt^ 
transition is determined by the well known constant gp^mr- This contribution is the dominant 
one for B — > tttt decays due to a large probability of i? — ^ pp-transition. Let us note that 
in the limit Mb ~^ 00 the ratio Br{B(i -^ pp)/Br{Bd -^ tttt) grows as M^, that is why FSI 
phase ^2 (pp) (and 5q (pp)) diminishes only as I/Mb- On the contrary tttt intermediate state 
plays a minor role in B ^ pp-decays. 

In description of tttt elastic scattering amplitudes in Eq.(12) the contributions of P, f and 
p Regge-poles were taken from ref.[16]. Finally vrai intermediate state should be taken into 
account. The large branching ratio of -B^ -^ Tr^a^^-decay ( Br(i?rf -^ n'^af) = (40 ±4)* 10^^) 
is partially compensated by the small pnai coupling constant (it is 1/3 of pTrvr one). As a 
result the contribution of irai intermediate state (which transforms into tttt by p-trajectory 
exchange in the t-channel) to FSI phases equals approximately that part of tttt intermediate 
state contribution which is due to p-trajectory exchange. Assuming that the sign of the Trai 
intermediate state contribution to phases is the same as that of the elastic channel and taking 
into account that the loop corrections to B ^ titi decay amplitudes lead to the diminishing 
of the (real) tree amplitudes by ^ 30% we obtain: 

S^ = 30° , 6^ = -10° , S^-S; = 40° . (14) 

The accuracy of this prediction is about 15°. 

For pp final state the analogous difference is about three times smaller, 6q — 62 ^ 15°. 
Thus the proposed model for FSI enables us to explain the B -^ tttt/B -^ pp puzzle. 

4 Direct CPV in B ^ TTTT-decays and phases of the 
penguin contribution 

It follows from Eq. (2) that the direct CP asymmetry in Ba{Bfi) -^ n^ii^ decay has the 
following expression in terms of quantities Aq,A2,P and phases: 

C+_ = -1= sinQ;[v/2Ao sin((5o - So - 6p) + A2 sm{62 - So - 6p)]/ 
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Thus the direct CP-violation parameter is proportional to the modulus of the penguin 
amplitude and is sensitive to the difference of the strong phases of Aq, A2 and penguin 
amplitudes. So far we have discussed the phases of the amplitudes Aq, A2. The penguin 
diagram contains a c-quark loop and has a nonzero phase even in the QCD perturbation 
theory. It was estimated in ref.[l] and is about 10°. Note that in PQCD it has a positive 
sign. 

Let us estimate the phase of the penguin amplitude 5p considering the charmed mesons 
intermediate states: B -^ DD, D*D, DD*, D*D* -^ tttt. In Regge model all these amplitudes 
are described at high energies by the exchanges of Z^*(D2)"^^3.jectories. An intercept of 
these exchange-degenerate trajectories can be obtained using the method of [17] or from 
the masses of D*(2007) - 1~ and D2(2460) - 2+ resonances, assuming linearity of these 
Regge-trajectories. Both methods give q;d*(0) = —0.8 -. — 1 and the slope a^, ~ 0.5GeV^^. 

The amplitude of D^D^ -^ n^ii^ reaction in the Regge model proposed in paper [18] 
can be written in the following form: 

TnD^.As,t) = -|e— Wr(l - az>*(t))(s/s,,)--W , (17) 

where r(a;) is the gamma function. 

The t-dependence of Regge-residues is chosen in accordance with the dual models and is 
tested for light (u,d,s) quarks. According to [18] Sd ~ 2.2 GeV^. 

Note that the sign of the amplitude is fixed by the unitarity in the t-channel (close to 
the Z)*-resonance). The constant Qq is determined by the width of the D* -^ Dir decay: 
(y'o/(167r) = 6.6. Using eq.(9) and the branching ratio Br{B -^ DD) ^ 2 ■ 10^^ we obtain 
the imaginary part of P and comparing it with the contribution of P in P ^ tt+tt^ decay 
probability we get 6p ~ —3.5°. The sign of 6p is negative - opposite to the positive sign 
which was obtained in perturbation theory. Since DD-decaj channel constitutes only f» 10% 
of all two-body charm-anticharm decays of Prf-meson, taking these channels into account we 
easily get 

6p 10° , (18) 

which may be very important for the interpretation of the experimental data on direct CP 
asymmetry. 

It was shown in ref.[2] that assuming that the phases satisfy the conditions: Sq — S2 = 
37°, 62 < and 6p > 0, it is possible to obtain the following inequality 

C+_ > -0.18 . (19) 

It is worthwhile to compare the obtained numbers with the value of C^ which follows 

from the asymmetry Acp{K~^7r~) ii d ^^ s symmetry is supposed [19]: 
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Experimental results obtained by Belle [20] and BABAR [21] are contradictory: 

C:f!_"" = -0.55(0.09) ,C:[^f^^^ = -0.21(0.09), (21) 

with Belle number being far below (19). For a non-perturbative phase of the penguin contri- 
bution (18) the value of the theoretical prediction for C+_ can be made substantially smaller 
and closer to the Belle result. 

For direct CP asymmetry in Bd{Bd) -^ tt^tt^ decay from (30) we readily obtain: 

Coo = ~y o^ sin a[Ao sin(5o - Sq - 6p) - V2A2 sm{62 - Sq - Sp)]/ 
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Coo ~ -1.06[0.8sin((5o - 5o - bp) - 1.4sin(52 -l^- bp)\ ^ -0.6 . (23) 

This unusually large direct CPV (measured by |Coo|) is intriguing task for future mea- 
surements since the present experimental error is too big: 

Co^-P- = -0.48(0.32) . (24) 

Another CPV asymmetry measured in B(i{Bct) -^ nn decays S^ is sensitive to the 

unitarity triangle angle a. Let us first neglect the penguin contribution. Then from the 
experimental value 5'^^'^'" = —0.62 ± 0.09 [20, 21] we get: 

sin 2a^ = 5+_ , (25) 

a^ = 109° ± 3° . (26) 

The penguin shifts the value of a. The accurate formula looks like: 

S^- = [sin2a(|±^±^cos(5o-52))- 

A P ~ [2 ~ 

— — -^ sin a cos(52 — 60 — ^p) — \ i^AqP sin a cos(5o — 60 — ^p)]/ 
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The numerical values of a from different B-decays will be given in the next Section. 



5 Analysis of Bci{Bci) -^ p^n^ decays 

The time dependence of these decay probabihties are given by the following formula [4]: 

dAt ^ ^^^ 

X [1 - q{Cp^ ± ACp^) cos(Amt) + q{Sp^ ± ASp^) sin(Amt)] , (28) 

where q = —1 corresponds to the decay of a particle which was B^ aX t = 0, while q = 1 
corresponds to the decay of a particle which was Bd at t = 0. According to [4]: 

|/1+-|2 _ M-+|2 I M+-|2 _ I /1-+|2 

where A^^ are the amplitudes of B^, -^ p'^vr^ decays, while A^^ are the amplitudes of 
Bd -^ p^TT^ decays. Introducing the ratios of the decay amplitudes: 

where q/p = e^^*'^ comes from B^ — Ba mixing and f3 is the angle of the unitarity triangle, 
we obtain the expressions for the remaining parameters entering Eq. (28): 
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where Cp.,^ and 5*^^ (as well as ^cp) ^^^ CP-odd observables, while ACp^^ and ASp.,^ are 
CP-even. The experimental data for the observables entering Eq. (28) accompanied by the 
averaged branching fraction are presented in Table 2 [3]. 

The decay amplitudes A^^ are described by the tree and penguin Feynman diagrams 
shown in Fig.l. The analogous diagrams describe amplitudes A^^. The corresponding 
formulas for the amplitudes look like: 



A-+ = Aie-''' + Pie^(^+''i) 
A-+ = A2e'^ + P^e-'^f^-^^^ 
A+- = A2e-'^ + P2e*('^+^2) 
A+- = Aie'^ + P-^e-'^f^-^'^ 

A1/A2 = ai/a2e'^ 
Pi/P2=pi/P2e'' , (32) 

where 7 and j3 are the angles of the unitarity triangle, while 61 and 62 are the difference of 
FSI strong phases between penguin and tree amplitudes (for penguin amplitudes we use the 
so-called t-convention, subtracting charm quark contribution to penguin amplitudes). 

All in all we have seven parameters in Eq.(32) specific for pvr final states (oi, a2,Pi,P2,Si, ^2 
and S) plus UT angle a = vr — /3 — 7, while the number of the experimental observables in 
Table 1 is six. To go further we should involve additional theoretical information in order to 
reduce the number of parameters. If we find the values of pi and p2 even with considerable 



uncertainties it will be very helpful for determination of UT angle a, since penguin ampli- 
tudes shift a by small amount proportional to Pi/ai, and even large uncertainty in this shift 
leads to few degrees (theoretical) uncertainty in a (see below). 

The most straightforward way is to calculate the matrix elements of the corresponding 
weak interactions Lagrangian with the help of factorization, as it was done in [22]. However 
it was shown above that there are substantial deviations from factorization in i? -^ vrvr 
decays. In particular from the experimental data on direct CP-asymmetry in Bfi{B(i) — > 
TT+TT^ decays we know that the factorization strongly underestimates the contribution of a 
penguin diagram to the decay amplitude [1, 2]. Another approach is to extract the penguin 
amplitudes from the branching ratios of the i?" -^ K^*Ti^ and B~ — > K^p'^ decays in which 
the penguin dominates with the help oi s ^^ d quark interchange symmetry, analogously to 
what was done for penguins in _B — i> tttt [23] and B ^ pp [24] decays. 

Feynman diagrams responsible for these decays are shown in Fig. 2. Comparing Fig. 2 
with Fig. 1 (b) we readily get the following relations: 
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1.071 
from which the values of pi and p2 follow: 

pI = 0.40(4) ■ 10~" , pI = 0.30(6) ■ 10"" , (35) 

where here and below we neglect the common factor IGirmBTB^, to which squares of am- 
plitudes are proportional. The remaining 8 — 2 = 6 parameters entering Eq.(32) we will 
determine from six experimental numbers presented in Table 2. 

From Eq.(32) we get the following relation for the averaged branching ratio of Bd{Bd) 
decays to p^vr^: 

ai+a| ^ Pi+P| ^ 23.1(2.7) ■ 10"" , (36) 

where the penguin-tree interference terms are omitted (being proportional to cos(7r— /3— 7) = 
cos a they are very small since UT is almost rectangular, a ~ vr/2). 
To determine the values of Oj the equation for ACp^^ is helpful: 

and from (35) - (37) and the experimental value for ACp^^ from Table 1 we get: 

a? = 31(3)-10~" , a^ = 14(3) ■ 10"" . (3^ 
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Now from the equations for Cf„ and Aq^ using the experimental data from Table 1 we 
are able to extract FSI phases 5i and 82'- 

2piai sin 5i + 2^202 sin (52 , al-al . , 

CpT, = 2 I 2 ^ / 2 I 2\2 [^^2^2 sm 62 - 2piai sm di\ , 

p^ 2piai sin 6i - 2^202 sin 62 , . 

Acp - 9 2 ' l-^yj 

sin(5i = -0.55(30) , sin^s = 0.51(40) , (40) 

and we see that large experimental errors of Cp^r and A^p do not allow the accurate deter- 
mination of the values of FSI phases. 

From the equations for S and AS* we will determine the values of a and 6: 

S,^ + ASp^= (41) 

aia2 sin(2Q; — 6) — pia2 cos(5i — 6) — ^2^1 cos(52 — S) + 2pia2 sin 61 sin 6 
of + al + 2piai sin 61 — 2^202 sin 62 

S^ - AS^ = (42) 

aia2 sin(2a + 6) — ^2^1 cos(52 + S) — pia2 cos(5i + 6) — 2p2ai sin 62 sin 6 
of + al + 2p202 sin 62 — 2piai sin 61 

where in (small) terms proportional to Pi we have substituted a = 7r/2. 

Substituting the numerical values for the parameters in the denominators we get: 

[(6±4)5p^-(45±4)AS'p^]10"^ = 2aia2sin5cos2a , (43) 



[(45 ± 4)5"^^ - (6 ± A)ASp^] 10"^ = 2aia2 sin 2a cos 5 - 

—2p2ai cos((5 — ^2) — 2pia2 cos(5 + 61) . (44) 

From the first equation we see that 6 equals zero or ir with ±5" accuracy. For UT angle 
a from the second equation neglecting the penguin contributions we obtain: 

aj; = 90° ± 3°(exp) , (45) 

while taking penguins into account we get: 

ftp^ = 84° ± 3°(exp) , (46) 

where 61 ~ —30° and 62 ~ 30° were used. 

Thus penguins shift a by 6° and even assuming 50% accuracy of c? ^^ s symmetry which 
was used to determine the numerical values of Pi allows us to determine ap.,^ with theoretical 
accuracy which equals the experimental one, originating from that in S'p^ and pointed out 
in (46): 

ap^ = 84° ± 3°(exp) ± 3°(theor) . (47) 
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The consideration of Bd{Bd) -^ tttt decays (see Eqs.(25) - (27)) leads to the following 
result: 

a^^ = 88" ± 4°(exp) ± 10°(theor) , (48) 

where a relatively large theoretical error is due to big (20°) shift of the tree level value of 
q;^^ by poorly known penguins and this time (unlike in [2]) we suppose 50% theoretical 
uncertainty in the value of penguin amplitude. 

In the case of Bd{Bd) -^ P^ P" decays penguin shifts the value of a by the same amount 
as is considered in this paper for Bd{Bd) -^ p^Ti^ decays, so the theoretical uncertainty is 
the same: 

ttpp = 87" ± 5°(exp) ± 3''(theor) , (49) 

while larger experimental uncertainty is due to that in Spp, 

S^p = -0.06 ± 0.18 , (50) 

which is twice as big as in Sp-,^. 

It is interesting to compare the numerical values (46), (48), (49) with the recent results 
of the fit of Unitarity Triangle [25, 26]: 

^CKMfitter _ ggo ^ go ^ (g^^, 

^UTfit _ g^o ^ go ^ (52) 

Large New Physics(NP) contribution to b ^ dg penguin could help to avoid large FSI 
phases since now the enhancement of direct CPV seen in A(^p will originate from closeness of 
tree level and penguin amplitudes. Also puzzle of large BrBd{Bd) -^ Tr^vr" can be resolved by 
NP contribution to b —^ dg penguin comparable with SM one recalculated from Bu -^ K'^ir'^ 
decay. The bound on such contribution comes from the coincidence within the errors of the 
values of a extracted from B ^' nn , pn and pp decays, where the penguin contributions are 
very different. ""^ 

These are strong arguments in favor of the measurements of the parameters oi B ^ tttt, 
pn and pp decays with better accuracy, which can be performed at LHCb and Super B 
factory. A search of NP manifestation by different values of UT angle a extracted from 
B ^> -n-n and B -^ np, pp decays is analogous to the one suggested in [27] through the 
difference of a extracted from the penguin polluted B ^ Tin decay and from UT analysis 
based on tree dominated observables Vcb and 7. 

At the end of this section let us note that the results (48) and (49) were obtained in the 
analysis based on isotopic invariance of strong interactions from the violation of which the 
additional uncertainty in a could follow [28]. Fortunately since in the absence of penguin 
amplitudes the relation S'^7r,pp = sin 2q;-^ is free from this type of uncertainty, it is manifested 
only as several percent correction to the shift of a induced by penguin which is negligible 
even for B ^ nn decays. 



^The same argument can be applied against large NP contributions to b —^ sg penguin: if the same NP 
does not enhance b ^ dg penguin the value of a from B -^ nn data will be closer to a^^ — 109° and disagree 



with that from a^rp and Upp. 



11 



6 Direct CPV in B ^ nK decays 

Recently Belle has published new results of the measurement of CP asymmetries in B(i{B(i) 
K+TT-{K-n+) and B+{B-) -^ K+TT^{K-n^) decays [29]: 



T{B, ^ K-n+) - T{Ba ^ K 



+ T,-) 



-^-('^'"'"' ^ r(B:^if-..) + r(B:^K..-) = -«-mms) , (53) 



T{B- -^ K-ttO)+T{B+ ^ K+tt^) 
In [29] the 4.5 standard deviations difference of these asymmetries was considered as a 



Acp{K+7r') = -A— ^^ ' ^ ^^^ -—L = 0.07(3)(1) . (54) 



paradox in the framework of the Standard Model (see also [30]) which it really were IF one 
neglected the color suppressed tree quark amplitude. Taking into account QCD penguin 
diagram and tree diagrams one easily gets the following relation between CP asymmetries 

[31]: 

Acp{K+n-) = AcpiK+n"^) + AcpiK'^n') , (55) 

where AcpiK^Tf^) is proportional to the color suppressed tree amplitude C. The exper- 



imental value of Acp{K tt ) has large uncertainty: 



Acp(i^%") = -0.14 ±0.11 , (56) 

however with the help oi d ^^ s interchange symmetry it can be related with CP asym- 
metry Coo of Bd{Bfi) -^ TT^TT^ decays: 



0_0\ 



AcpiK'n 



r B, -^ tt^tt'J + r B, -^ Ti^Ti' 



V.,..M 



us ^ ts 



Vt 



td 



"""C™ . (57) 



sma 



where the opposite signs in the definitions of Aqp and Coo s'-re compensated by a negative 
sign of Vts- The experimental uncertainty of Coo is also very large, that is why we use the 
above result (Eq.(23)) for numerical estimate: 

Coo ~ -0.6 . (58) 

Substituting (58) in (57) and (53), (54) and (57) in (55) we finally obtain: 

-0.094 ±0.02 = (0.07 ±0.03) + (-0.07 ±0.02) , (59) 

resolving in this way the paradox noted in [29] (the remaining ^ 2a difference can be 
safely attributed to statistical fluctuation). Concluding this Section let us remind that the 
absence of color suppression of the tree amplitude of Bd -^ tt'^tt*' decay is explained in 
Sections 2,3 by large FSI phases difference of tree amplitudes with isospin zero and two. 

7 Polarizations of vector mesons in B ^ l/l/-decays 

In this Section we consider Bd{Bd) decays into the pair of light (p, K*, if) vector mesons. The 
short distance contributions to vector meson production in B-decays lead to the dominance 
of the longitudinal polarization of the vector mesons. This is a general property valid in 
the large Mq- limit due to helicity conservation for vector currents and corrections should 
be ~ My/Mq. It is satisfied experimentally in B ^ P^P^ decays, where the contribution 
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of longitudinal polarization of p mesons is Jl = T^/F = 0.968 ± 0.023. Let us note that 
FSI are not important for these decays; for example there are no large strong interaction 
phases generated by rescattering. The manifestation of this statement is the absence of 
enhancement of the color suppressed amplitude which describes the decay into p°p°. 

On the other hand there are several B-decays to vector mesons, where the longitudinal 
polarizations give only about 50% of decay rates. For example: 

for B+ -^ K*^p+ Jl = 0.48 ± 0.08, B^ -^ K*^p^ Jl = 0.57 ± 0.12, B+ -^ (j)K*+ Jl = 
0.50 ± 0.07, Bd -^ (pK*° /l = 0.491 ± 0.032 [3]. 

This is a real puzzle IF only short distance dynamics for these decays is invoked. We would 
like to argue that strong rescattering related to large distance dynamics may be responsible 
for the observed polarizations pattern. First let us note that in all the decays, where Jl ~ 
50%, the penguin diagrams give dominant contribution. In this case a large contribution to 
the matrix elements of the decays comes from DDs{D*Ds, DD*, ..) intermediate states, which 
have large branching ratios. In Section 4 we analyzed DD intermediate state contribution 
to a strong phase of the penguin amplitude for BdiBd) -^ tttt decays. It was argued that 
the phase of the order of 10° can be generated by the charmed mesons intermediate state. ^ 
FSI is important due to large branching ratio Br(i?rf — > DD) ^ 2 ■ 10~^ in comparison with 
the penguin contribution to decays to two pions [2]: Br(i? -^ 7r7r)p = 0.6 • 10~^ extracted 
with a help oi d ^^ s symmetry. In the case of b ^ sg penguin dominated decays the 
intermediate state contains DDg pair, and we should compare Br{Bd -^ DDg) ~ 10~^ with 
Br{B -^ K*^p'^) K, 10'^. That is why the relative contribution of DDs states are ~ 2 
times larger than for B -^ nn penguins. The amplitude of the binary reaction DDg -^ VV 
at high energies is dominated by the exchange of D*-regge trajectory and according to 
general rules for spin-structure of regge vertices (see for example [32]) the final vector mesons 
are produced at high energies transversely polarized. Thus we expect a large fraction of 
transverse polarization of vector mesons in these decays. The value of /x, is sensitive to 
intercept of D*-trajectory [33]. If the penguin contribution in the decays indicated above is 
dominant in the SU(3) limit we have: 

Br{(j)K*^) = Br{K*^p'^) = Br{(t)K*+) = 2Br{K*^p^) (60) 

and fi in all these decays should be the same. These predictions agree with experimental 
data [3]. 

8 Puzzle of charm- ant icharm baryons production 

Large probability of B-decay to AcSc has been observed recently: Br{B^ -^ ^f^^c ~ 10~^) 
[3]. It is surprisingly large compared to the branching of B-decay to A+p = (2.19 ±0.8)10"^. 
From PQCD point of view both processes are described by similar diagrams with the sub- 
stitution of ud (for p) by cs (for Sc) and phase space arguments even favor p-production. 

On the other hand from the soft rescatterings point of view the large probabilities of 
DDs{D*Ds,DD*, ...) intermediate states, considered in the previous section, can play an 
important role in B^ -^ A^S°-decays. For A+p final states the corresponding two-meson 
intermediate states have smaller branchings and, what is even more important, have different 
kinematics. For DDg, ... intermediate states the momentum of these heavy states is not large 
{p ^ 1.8 GeV) in B rest frame and all light quarks {u, d, d, s) are slow in this frame. The final 



^Let us emphasize that this considerable contribution to the subdoniinant penguin aniphtudc is not 
important for the B^(B^) -^ tttt widths, where the tree diagram dominates. 
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A+S° are also rather slow in the B-rest frame and thus all quarks have large projections to 
the wave functions of the final baryons. On the contrary for 7rD,pD, .. intermediate states 
in A+p-decays momenta of u, d- quarks in light mesons are large and the projections to 
the wave functions of final baryons have extra smallness. The resulting suppression can be 
estimated in regge-model of ref.[18] with the nucleon trajectory exchange in the t-channel 
and is ~ 10~^ in accordance with experimental observation. 

9 Conclusions 

FSI play an important role in two-body hadronic decays of heavy mesons. Theoretical 
estimates with account of the lowest intermediate states give a satisfactory agreement with 
the experiment and provide the explanation of several puzzles observed in B decays. 
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Table 1 



Mode 


Br(10-6) 


Mode 


Br(10-6) 


Bd -^ 71 + 71- 


5.2 ±0.2 


Bd -^ P^P~ 


24.2 ±3.2 


Bd ^ ttOttO 


1.3 ±0.2 


Bd ^ pV° 


0.68 ±0.27 


Bu -^ tt+ttO 


5.7 ±0.4 


Bu ^ P+P° 


18.2 ±3.0 



C-averaged branching ratios of i? ^ tttt and B ^ pp decays. 
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Table 2 



BrB,{Ba) ^ 

^p±7rT 




C 


AC,. 


^piT 


A5,. 


(23.1 ±2.7)10^6 


-0.13 
±0.04 


0.01 
±0.07 


0.37 
±0.08 


0.01 
±0.09 


-0.04 
±0.10 



The experimental values of observables which describe Bd{Bd) —* p n^ decays 
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FIGURES 
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Figure 1: Tree and penguin diagrams for B-decays to pvr- mesons. Bd ^ p vr"*" decay is 
described by the amplitudes Ai and Pi, while Bd -^ p^vr^ decay - by the amplitudes A2 and 

P2- 




^ 



o* ^o 



^" r 






Figure 2: B decays in which the penguin diagram dominates. 
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